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ABSTRACT 

Most galaxies follow well-defined scaling relations of metallicity (O/H), star formation 
rate (SFR), and stellar mass (M s t ar ). However, low-metallicity starbursts, rare in the 
Local Universe but more common at high redshift, deviate significantly from these 
scaling relations. On the "main sequence" of star formation, these galaxies have high 
SFR for a given Mst ar ; and on the mass-metallicity relation, they have excess M s t ar for 
their low metallicity. In this paper, we characterize O/H, M star , and SFR for these de- 
viant "low-metallicity starbursts", selected from a sample of ~1100 galaxies, spanning 
almost two orders of magnitude in metal abundance, a factor of ~ 10 6 in SFR, and of 
~ 10 5 in stellar mass. Our sample includes quiescent star-forming galaxies and blue 
compact dwarfs at redshift 0, luminous compact galaxies at redshift 0.3, and Lyman 
Break galaxies at redshifts 1-3.4. Applying a Principal Component Analysis (PCA) 
to the galaxies in our sample with M star ^ 3 x 1O 1O M0 gives a Fundamental Plane 
(FP) of scaling relations; SFR and stellar mass define the plane itself, and O/H its 
thickness. The dispersion for our sample in the edge-on view of the plane is 0.17 dex, 
independently of redshift and including the metal-poor starbursts. The same FP is fol- 
lowed by 55 100 galaxies selected from the Sloan Digital Sky Survey, with a dispersion 
of 0.06 dex. In a companion paper, we develop multi-phase chemical evolution models 
that successfully predict the observed scaling relations and the FP; the deviations 
from the main scaling relations are caused by a different (starburst or "active" ) mode 
of star formation. These scaling relations do not truly evolve, but rather are defined 
by the different galaxy populations dominant at different cosmological epochs. 

Key words: galaxies: abundances - galaxies: dwarf - galaxies: evolution - galaxies: 
high-redshift - galaxies: starburst - galaxies: star formation 



1 INTRODUCTION 

Our understanding of how galaxies assemble their mass over 
cosmic time has advanced considerably in the last decade. 
There have emerged observationally-defined scaling relations 
that link metallicity, stellar mass, and star formation rates 
(SFRs) in galaxies for lookback times of up to ~10Gyr. 
The rate at which galaxies form stars is correlated with 
their stellar mass, up to > 1O 11 M where the trend flattens. 
This trend, known as the "main sequence of star formation" 
(SFMS) has roughly the same slope up to redshifts z ~3, but 
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an increasing normalization such that higher redshift galax- 
ies produce more stars for their mass than low-redshift ones 
jNoeske et al l 12007): " 



l2007t iBauer et all (20111 ). Thus 



Sali m et all 120071 : | Schi minovich et al.l 
high-redshift populations 



tend to have higher specific SFRs (SFR divided by stel- 
lar mass, sSFR) than nearby galaxies, although there is still 
some debate about whether at all redshifts higher sSFRs are 
observed in less m assive galaxies than in their more massive 
coun t erparts (e.g.. IZheng et al.1 120071 : iPerez-Gonzalez et al.l 
|200S| ; iRodighiero et al.ll2010D . 

Metallicity and stellar mass are also related (the 
mass-metallicity rela tion, MZR), but through a somewhat 
tighter correlation (|Tremonti et al.1 120041 : ISavaglio et al.1 



L. K. Hunt et al. 



20051; lErb et alJlioO^ ; iMaiolino et"ai1l2008l ; iMannucci et all 



20091 ). A secondary dependence of the MZR on SFR, 
dubbed the "Fundamental Metallicity Relation" (FMR), re- 
duces the scatter cons i derably for local galaxy populations 
l|Mannucci et alj|20ld ). IMannucci et all (|201lh extended the 
FMR to lower stellar masses with Gamma-Ray-Burst host 
galaxies, but it is still unable to encompass the properties 
of galaxies at z > 3, and as we show below, also in nearby 
metal-poor starbursts. 

Although the SFMS and the MZR (and FMR) are good 
at predicting the links among star formation, metallicity, 
and stellar mass for the bulk of galaxies up to redshifts of ~2, 
some subsets of galaxy populations show clear departures 
from the main trends. In particular, vigorous starbursts with 
high SFRs tend to lie above the SFMS at all redshifts and 
to have too much mass (or luminosity) for their metallic- 
ity (e.g..lHovos et al. 2005; R osario et a l. 2008; S alzer et al.l 



2009 ; Peeples et alll2009l : iRodighiero et al.ll201ll ; lAtek et al 



201 If ). Usually, these high-SFR outliers are attributed to the 



merger mode of star formation; galaxy-wide SF occurring 
over long dynamical times is thought to define the SF Main 
Sequence, while merger-induced bursts of star formation re- 
sult in outliers from the main trend (e.g.. IRodighiero et al.l 
l201ll ; lLee et~"ai1l2004l ). 

In a companion paper l|Magrini et alll2012ft . we pro- 
pose that mergers may not be directly responsible for pro- 
ducing the starbursts that deviate from the general SFR- 
mass-metallicity relations, but rather that initial physical 
conditions, independently of their origin, may be driving 
the deviations. We develop different classes of models based 
on distinct sets of initial conditions, including mass and ini- 
tial density of the galaxy, and size and density of the star- 
forming regions. These would correspond to different modes 
of star formation, an "active" or starburst mode, and a more 
"passive", quiescent one, with which we successfully model 
the scaling relations of <1100 galaxies. Here, we present the 
sample of galaxies, and describe the observed scaling rela- 
tions of metallicity, stellar mass, and star formation rate for 
the sample. 



1.1 Our approach 

In this paper, we want to investigate the empirical connec- 
tion between high-z starbursts and metal-poor galaxies with 
high SFRs in the Local Universe. Our aim is not to analyze 
the bulk of the star-forming galaxies representative of the 
Local Volume, as has already been done by several groups 



with the SPSS (e.g., 



emonti et al. 2004; Mannucci et al. 



2010| : lLara-Lopez et alJl201ol : lYates et al.ll2012l ; IZahid et all 



20121 ). Rather, our intention is to identify classes of galaxies 
which deviate from the SFMS and the MZR in the same way, 
placing particular emphasis on potential local analogues of 
high-redshift galaxy populations. The idea is that the out- 
liers of well-defined scaling relations may reveal insights that 
cannot be gained otherwise. Hence, we require a sample that 
spans a wide range in parameter space, in order to ade- 
quately characterize trends with metallicity, SFR, and stel- 
lar mass over a large redshift interval. 

As we show below, some low-metallicity Blue Com- 
pact Dwarf galaxies (BCDs) have high SFRs, particularly 
high sSFRs, and occupy the same locus of the MZR and 
the SFMS as z = 3 Lyman Break Galaxies (LBGs). Rel- 



ative to the main trend, these objects have an excess of 
>100 in stellar mass at a given metallicity, and a factor of 
~ 10 excess SFR at a given mass. The same is true for 
some of the "Green Peas" and Lu minous Compact Galaxies 
(LC Gs) at z ~ 0. 1 -0.4 selected bv lCardamone et all (|2009h 
and llzotov et al.l (|201ll ) from the Sloan Digital Sky Survey 
(SDSS). It seems evident that SFR (or SFR surface density 
Ssfr, which is more difficult to measure) is governed by 
similar physical conditions at all redshifts, and that which 
galaxy populations are observed at a given redshift depends 
on how they are selected. Because extreme starbursts are 
rare locally, but more common at high redshift, selection ef- 
fectsQ may be paramount in defining observed trends among 
stellar mass, metallicity, and SFR. 

The paper is organized as follows: in Sect. [2] we present 
several samples of dwarf galaxies in the Local Universe, two 
samples of intermediate- redshift galaxies (the "Green Peas" 
and LCGs), and several samples of LBGs at z ~ 1 — 3. We 
also explain how we derive stellar masses from 4.5 [im ob- 
servations for the nearby galaxy samples. Sect. [3] identifies 
a class of "low-metallicity starbursts" on the basis of devi- 
ations from the main scaling relations observed at low and 
high redshift. A "Fundamental Plane" of metallicity, stel- 
lar mass, and star formation rate independent of redshift is 
introduced in Sect. [4] and its characteristics are compared 
with common perceptions of the SFMS and the MZR. In 
Sect. [5] we discuss possible similarities between nearby star- 
bursts and starbursts at high redshift, and possible ramifi- 
cations for our current conception of galaxy evolution. 



2 OBSERVATIONAL SAMPLES 

Because of the need to compare SFR, metal abundance [as 
defined by the nebular oxygen abundance, 12+ log(0/H)], 
and stellar mass, M s t a r, we have selected only samples of 
galaxies which either have these quantities already avail- 
able in the literature, or for which we could derive them 
from published data. Five sets of galaxies in the Local 
Universe met these criteria : the nea r by dw arf irregular 
galaxies (dlrr) published by iLee et all (|2006l . here dubbed 
"dlrr"), the 11 Mpc d istance-limited sample of nearby galax- 
ies (T1H UGS. LVL: iKennicutt et al] 120081 ). the starburst 
sample of lEngelbracht et al.l d2008l), and the tw o BCD sam- 
ples s tudied by Fumagall i et al] (|201Ch . and iHunt et aD 
(|2010h . We added to these seven galaxy samples at higher 
redshifts: the "Green Pea" compact gal axies identified by 
the G alaxy Zoo team at z ~ 0.1 — 0.3 l|Cardamone et aD 
2009h. the LCG sample at z ~ 0.1 - 0.6 selected b y 
Izotov et all d201lt>. LBGs a t z ~ 1 dShaplev et all l2005ah , 
z ~ 2 llShaplev et all l2004l; lErb et al l 120061 ). and 2 ~ 3 
l|Maiolino et all |200S| ; IMannucci et all l2009h When stellar 
masses were not available in the literature, we determined 



1 We use the term "selection effect" to mean "the ten- 
dency for a conclusion based on observations to be influenced 
by the method used to select objects for observation" (see 
http://earthguide.ucsd. edu/virtualmuseum/ 'Glossary- Astro/ 
glosses- z.shtml). In this sense, how a galaxy (or galaxy popula- 
tion) is selected may be important for defining the slope and nor- 
malization of scaling relations, both locally and at high redshift. 
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stellar masses from IRAC photometry as discussed in Sect. 
12.61 Altogether, we consider in the analysis 1070 galaxies 
from z ~ to z ~ 3. Except for the high-redshift samples 
JShaplev et alj|2004 l2005al ; lErb et al.|[2006l ; iMaiolino et ail 
120081; Mannucci et al. 200£ ) and some of the LVL galaxies 
|Marble et al.ll2010l ). metallicities for all galaxies are deter- 
mined by the "dir ect" (electron temperature) method (e.g., 
Ilzotov et al. 1 12007ft since some of the samples are defined by 
requiring detections of [Om]A4363. We adopt the central 
abundance of each galaxy, and thus abundance gradients, if 
present, are not taken int o account . All s tellar masses and 
SFRs are reported to the IChabrierl |2003T ) scale (see below 
for specific instances). Each sample is described in some de- 
tail in the following. 



2.1 The dlrr sample 

I Lee et ail (|2006t ) studied 25 dwarf galaxies within 5Mpc 
distance that were observed at 4.5 urn with Spiteer/IRAC. 
The advantage of using 4.5 vim to measure stellar masses 
is that the stellar mass-to-light (M/L) ratio does not de- 
pend stro ngly on ag e or m etallicity, and dust extinction is 
negligible. I Lee et al.l ([200(a ) derived M/Ls for their sample 
using the 4.5 urn luminosity, and the observed B— [4.5] color, 
transforming the 4.5 urn (Vega) magnitude, [4.5], to K as- 
suming if— [4.5 ] = +0 .2. Sub-solar metallicity models from 
iBell fc de Jond i|200ll ) were then used to derive M/L(K) 
as a linear function of B — K; the solar absolute magni- 
tude at 4.5 urn was taken to be ~ +3.3, which implies that 
K — M color is roughly 0. All stellar masses are scaled to the 
IChabrierl (120031 ) Initial Mass Function (IMF). Oxy gen abun- 
dances for 21 galaxies are derived with the direct method 
when electron temper atures are availab le from [Oiii] A 4363 
measurements, or the iMcGaughl (|l99lh bright-line calibra- 
tion for three additional o bjects. 

SFRs are taken from IWoo et al.l (|2008l ) and iLee et all 
(2009), based on extinction-corrected Ha luminosities us- 
ing the iKennicuttl |l998) conversion factor. For 5 galaxies 
these were unavailable, so we adopted the SF R values in- 
ferred from the star-formation histories given bv lWeisz et al.l 
<|201ll) . 

A linear fit to the MZR for this nearby dwarf sam- 
ple gives a small scatter of 0.12 dex, over a wide range 
of stellar masses and metalli cities: ^10 6 to 10 9 Mq and 
7.35 (Leo A) to 8.3 in O/H (|Lee et al.|[200r3 ). The SDSS 
MZR has a comparable scatter, but does not extend either 
to low stellar masses (> 10 8 5 Mq) or to low metallicities 
(12+log(0/H)>8.3) (|Tremonti et al.ll2004T ). 



2.2 The 11HUGS/LVL sample 

iKennicutt et al.l l|2008l ) defined a distance-limited sample 
of 261 late-type galaxies (in the "primary sample") within 
UMpc, at high Galactic latitude (\b\ ^ 20°), and with 
Hubble type T > 0. They also imposed a magnitude limit, 
B <15. The UMpc Ultraviolet and Ha Survey (11HUGS) 
and the Spitzer Local Volume Legacy project (LVL) sub- 
sequently obtained Galaxy Evolution Explorer (GALEX) 
imaging, Ha photometry, and infrared (IR) imaging for the 
sample. Because the LVL is a fair sampling of the Local 
Universe, it is dominated by dwarf galaxies (see below). 



The stellar masses have been derived by us from 
the 4.5 urn lumin osities using the photometry given by 
Dale et ail (120091 ). We have adopted the formulation of 
Lee et al.l (|2006r ). but take the contamination by nebular and 
hot dust emission into account as described in Sect. 12.61 Stel- 



lar masse s have also been estim ated for the 11HUGS/LVL 
sample bv lBothwell et all l|2009h from the B—V color and B- 
band luminosities, using the formulation of Bell fc de Jond 
(|2001r) for the dependence of M/L on coloiQ. Their masses 
range from 10 7 to > 10 11 Mq, similar to our values. 

Oxygen abundances are available for 129 of the 
11HU GS/LVL galaxies, and are tabulated in iMarble et al.l 
(2010)e. The values are from either [Oiii] A4363, or em- 
pirical strong-line methods, but have been analyzed by 
IMarble et alTl|2010l) and show no systema tic deviations such 
as discussed by iKewlev fc Ellison! (|200ST ). The SFRs have 
been taken from Lee et al. (2009), using the values derived 
either from extinction-corrected Ha or UV luminosities, 
whichever are larger. We considered both SFR values and 
chose the largest one beca use Ha may not be a robust tracer 
of SFR at low SFRs (e.g., iMelena et aT]|2009l ; iBoselli et all 
120091 ; II ee et al.ll201ll ). and we wanted to ensure an aggres- 
sive correction for nebular emission (see Sect. l2.6ll . There are 
a total of 113 11HUGS/LVL galaxies with all three quanti- 
ties (M s t ar , O/H, and SFR) available. 

The 11HUGS/LVL sample is dominated by dwarf galax- 
ies; 90 (78%) of these 1 13 LVL galaxies have Hubble types 
T ^ 7 (|Lee et al.ll2009h . corresponding to Sd or later, and 
79 (68%) with T ^ 8, late r than but including Sdm. As dis- 
cussed bv lLee et al.l (|2009h . there are a few dwarf spheroidals 
with 12+log(0/H)^ 8.4; 3 of these remain in our LVL 
subset. SFRs generally range from 3 x 10 -4 Mq yr _1 to 
~ 7M yr- 1 ; NGC 784 is a metal-poor (12+ log(0/H)=7.9) 
outlier with 23M0yr _1 . Oxygen abundances run from 7.2 
(UGC 5340) to 9.3. There are some weak active galactic nu- 
clei in the LVL (e.g., NGC 3726, NGC 4826), but they are all 
metal rich and we have eliminated them from the analysis. 



2.3 The BCDs 

The dwarf galaxies most likely to be metal poor with 
high SFRs are the BCDs. These galaxies were defined 
originally on the basis of blue colors and compact mor- 
phology dThuan fc Martini Il98lh . but subsequently "com- 
pact" was superseded by high (stellar) surface density 
l|Gil de Paz et alj 120031 ). The parent samples from which 
most BCDs are drawn are generally objective prism surveys 
(e.g., Markarian, SBS, UM, Tololo, etc.) which tend to pref- 
erentially select nearby galaxies with high-equivalent-width 
emission lines. 

Our BCD sample is the union of the starbursts studied 
Engelbracht et al. (I2008T ). a nd the galax i es inv estigated 



Fumagalli et alj (|2010h and iHunt etail (|201Cj ). As de- 



by 
by 

scribed below, there are a total of 89 BCDs in the combined 
samples (excluding 5 galaxies already in the LVL). Not all 
the galaxies in the Engelbracht sample are BCDs or dwarf 
galaxies; it was originally designed as a "starburst sample" 



2 Because B — V observations were only available for roughly 
half the sample, they assumed a standard B — V color for the 
remaining galaxies. 
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with a wide range of metallicities. Hence, it adds high-mass 
star-forming galaxies to the mix. 

Included here are also the two prototypical "ac- 
tive" and "passive" B CDs, SBS 0335-052 and IZwl8 
|Hirashita fc Hu nt 2004). Although these two galaxies have 
very similar metallicities, 12+ log(0/H)~7.2, the ir SFRs dif- 
fer by a fac tor of 10 (1 M© yr" 1 vs. 0.1 M© yr" 



Hunt et al 



l200ll. [2005). We use the masses given by iFumagalli et al 
1 20101 ). which agree well with the masses we obtain here 
from the 4.5 \xm luminosity (see Sect. I2.6")l . 

The distinction between a BCD and an "active" dlrr 
is not clear cut, and may partially be a question of dis- 
tance and the se nsitivity with whic h extended emission 
can be observed jTolstov et alj|2009h. I n fact , not all the 
dwarf galaxies in the lEngelbracht' "et all (|2008h sample are 
bona-fide (high surface brightness, compact) BCDs. There 
are two duplications with the LVL sample (UGC4483 and 
UGCA292); in both cases, we use the LVL determina- 
tion of O/H, SFR, and stellar masses. On the other hand, 
some LVL galaxies are considered BCDs, e.g., NGC5253, 
NGC 1140. There are also some d uplications among the 
BCDs; Mrk 209 appears i n both the IFumagalli et alj l|201Ch 
and the lHunt etal] l|201Ch samples and II Zw40, NGC 2537, 
NGC 5253, UM462, UM 462 are i n both lEngelbracht et all 
J2008I) and IFumagalli etal] (|2010h . In all cases, we use the 
Fumagalli et al. l|20ld ) values for O/H and SFR. 

For the BCD samples, as for the 11HUGS/LVL sample, 
we derive stellar masses from the 4.5 \im luminosity, tak- 
ing nebular and hot-dust contamination into account (see 
Sect. 1 2. 6[) . Oxygen abundances a re taken from | lzotoy et al.l 
(2007), or from the compilation bv lEngelbracht et alT (2008). 
Virtually all of these are derived using [Oin] A4363, so should 
be re latively robust in the fa ce of systematic uncertainties 
fe.g.. lKewlev fc Ellison! [20081 ). 

SFRs were calculated from extinction-corrected Ha lu- 
minosities, or alternatively, when the former were unavail- 
able, from IR luminosities Lir a ccording to the formula- 
tion gi ven by Draine fc Lil lj2007h and the conversion fac- 



tor of Kennicutt |l998h . The IR photometry came from 
lEngelbracht et all (§008) or from Hunt et al. (2012, m prepa- 
ration). When both Ha and IR photometry were available, 
we took the largest value of SFR; in most cases (but not 
all), the Ha-derived SFRs are largest. 

There are several weak active galactic nuclei (AGN) 
with super-solar abundances in the Engelbracht et al. sam- 
ple, and we eliminated them from further consideration. The 
stellar masses would be uncertain because of possible hot 
dust from the AGN, and the abundances may also be incor- 
rect because of AGN contamination of the line emission. 

There are 89 galaxies in the combined BCD sample. The 
lowest abundance in this sample is 12+ log(0/H) = 7.1 for 
SBS 0335— 052W. In the metallicity range we are consider- 
ing, stellar masses (see below) range from 2x 10 5 Mq to 1.6 x 
10 n M Q , and SFRs range from 10" 4 M Q yr -1 (DDO 187) 
to 74M yr _1 (SHOC391, roughly 3 times higher than the 
largest LVL SFR). 

2.4 The Green Peas and the LCGs 

Recently, ICardamone et aL (2009) identified an unusual 
class of galaxies noted for their compact, green appearance 
on the SDSS composite images, the so-called "green pea 



galaxies" . These galaxies appear green because of their very 
strong [Om]A5007 optical emission line redshifted into the 
SDSS r band. The Green Peas have redshifts of 0.11 < z < 
0.36, and are low mass (10 8 5 - 10 10 M©) with high SFRs 
(~ lOMpyr" 1 ). T heir me tallicities were reex amined by 
lAmorfn et all (|2010h and by llzotov et all (|201ll ) and found 
to be somewhat metal poor, 12+log(0/H)^ 8.0. They are 
offset by ~0.3 dex in the MZR from the galaxies in the Local 
Universe (e.g.. iTremonti et al]|200"3 ). and are thought to be 
local analogues of more distant ultravio let-luminous galaxies 
such as LBGs (|Cardamone et alj|2009h . 



llzotov et al.l (l201ll ) selected a larger sample of sim- 
ilar galaxies, the LCGs, by introducing a threshold in 
(extinction-corrected) H/3 luminosity (3 x 10 40 ergs _1 ), re- 
quiring a high H/3 equivalent width (>50A), and including 
only those objects with a well-detected [Oiii]A4363 emission 
line in order to obtain ac curate abundances using the direct 
method (in contrast to ICardamone et all [2009, who used 
the less accurate empirical strong-line method). Moreover, 
a compact appearance was required on the SDSS images, 
and galaxies with evidence of AGN spectral features were 
excluded. These criteria resulted in a LCG sample of 803 
galaxies, with redshifts ranging from 0.1 to 0.6 (although 
most of them have 0.1 < z < 0.3), and metallicities from 
12+log(0/H) from 7.5 to 8.4, with a median value of ~8.1. 
SFRs are derived from extinction-corre cted Ha luminosities 
and range from 0.7 to ~ 6OM yr _1 . llzotov et all (|201lh 
compared these values with SFRs inferred from FUV lumi- 
nosities from GALEX, and found good agreement. 

The LCGs and the Green Peas provide excellent exam- 
ples of low-metallicity starbursts. They are at sufficiently low 
redshifts (2 ~ 0.1 — 0.3) that they should not present sig- 
nificant evolutionary differences relative to the less distant 
metal-poor starbursts described in the previous sections. On 
the other hand, their properties are sufficiently extreme in 
terms of low oxygen abundance and high SFR that they 
can be good local counterparts for hig h-redshift LBGs, a l- 
though they are much fainter (~ 3 mag, llzotov et al.ll201ll ). 
and less massive (see below). The SFR maximum value, 
~ 60 Mq yr -1 of the LCGs and Green Peas is virtually iden- 
tical to the maximum values (~ 60—66 M© yr" 1 ) in the LBG 
z ~ 2 and z ^ 3 sampl es described bv lSteidel et al.1 (|2004T ) 
and lSteidel et aL (2003). However, the minimum SFR in the 
high-redshift LBG samples is ~ 3 — 6Mq yr -1 , while in the 
LCGs is lower (by definition, because of the threshold in Ha 
luminosity), 0.7M©yr _1 . 



Here we include in our analysis the llzotov et al.1 (|201 lh 
LCG sample of 803 ga laxies, together with the subset of 
Cardamone et aL] (I2009T ) Green Pea galaxies (66) for which 



Izotov et al.l ( 201 ll ) remeas ured the mas s es and metal abun- 



dances (see their Table 2). llzotov et al.l l|201ll ) derived stel- 
lar masses by fitting the SDSS spectrum from 0.39 to 
0.92 u.m, and approximating the star-formation history with 
two episodes: a recent short burst and a previous continuous 
star-formation event r esponsible for t he evolved stars. They 
used a Salpeter IMF |Salpeteijll955l). and we have divided 
by a correction factor of 1.8 (see lLee et alj|2005 lErb et aL] 
120061; iPozzetti et al.ll2007r i to bring the masses and SFRs to 
the lChabrierl <|2003h scale. 
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2.5 The LBGs 

We use LBGs as a natural high-redshift galaxy popula- 
tion to compare with the local metal-poor starbursts. 12 
L BGs in the DEEP2 sa mple at z ~ 1 were characterized 
bv lShaplev et alj (|2005ah . The weakness of the [Om]A4363 
line precluded direct methods for measuring nebular oxygen 
abundance, so they used strong- line techniques (N2, 03N2). 
Stellar masses were calculated by fitting the BRIK S broad- 
band spectral energy di stributions (SEDs) with synthetic 
stellar populations from lBruzual fc Cha riot (200jl) and ex- 
ponentially declining SF histories. The y derived SFRs fr om 
Ha luminosities with the conversion o f lKennicutj (|l998l ). 

lErb et al.1 (|200rj ) studied a sample of LBGs at z ~ 2, 
and derived stellar masses by fitting the broadband SEDs us- 
ing U„GRJK , and, when available, Spit zer/IKAC bands to 
stand ard spectral synthesis models from lBruzual fc Chariot! 
(2003). SFRs were calculated from extinction-corrected Ha 
luminosities, and metallicities were established using the 
strong-line method because of the weakness of [On]A4363 
in these high-redshift objects. They arranged their sample 
in 6 bins of stellar mass, in order to ensu re rob ust trends 
with m ass and metallicitv. lErb et al.l |2006l ) use the lChabriej 
|2003l ) IMF, so we have applied no correction to the masses 
or the SF Rs. Another s a mple of 7 LBGs at z ~ 2 was 
studied bv lShaplev et al.l (I2004T). They fi t broadband SEDs 
(U„GRK S ) to iBruzual fe Charlotl (|2003h populations with 
a Salpeter IMF, so we have divided by a factor of 1.8 to 
correct them to the Chabrier IMF. As above, oxygen abun- 
dances were calculated from strong-line methods, and SFRs 

from Ha luminositi es. 

_ Maiolino et al.1 (|200Sl . AMAZE) and iMannucci et al.1 
(|2009l . LSD) studied 18 similarly selected LBGs at z ~ 3 but 
also required that at least two of the Spiteer/IRAC bands 
be present in order to obtain more reliable stellar masses. 
IMannucci et~ai1 (|2009l ) fit UGRJK + 3.6, 4.5 ^m broadband 
SED s with models from HYPERZ-MASS l|Pozzetti et al.1 
120071 ); HYPERZ-MASS considers a Chabrier IMF, so we 
have not introduced any corr ection factors for their ob- 
jects. In iMaiolino et all |2008l ). stellar masses are derived 
by fitting the broadband SED (with the same photome- 
try as above) with a s eries of standa r d spe ctral synthesis 
models as described by iFont ana et al l l|2006l) . They used a 
Salpeter IMF, so -as above- we have divided stellar masses 
and SFRs by a correction factor of 1.8 to be consistent with 
our adopted Chabrier IMF. Both papers estimate metallici- 
ties by comparing three independent strong-line calibrations 
of the direct method. 



2.6 Stellar masses from IRAC 4.5 /im 

For the dwarf galaxy samples in the Local Universe, we cal- 
culated the stellar masses from Spitzer/IKAC 4.5 vim lu- 
minosities , using a variation of the method described in 
I Lee et all (|2006l ). Near- infrared bands (A ^1 vim) in gen- 
eral, and IRAC bands in particular (because of their avail- 
ability and sensitivity) are particularly well suited for esti- 
mating stellar mass. At these wavelengths, the stel lar M/L 
ratio is quite insensit ive to age and metallicity |Gavazzil 
1 19931 ; [j un fe Iml 120081 . e.g.,), and dust extinction is negli- 
gible. However, nebular emission in metal-poor starbursts 
can be very significant, both in the continuum and in hy- 



drogen recombin ation and other lines (jReines et al.l l2010l : 
lAtek et al.ll201ll ). This can have a potentially strong im - 
pact on the 4.5 um photometry (|Smith fe Hancock] 12009). 
and thus on the stellar masses derived from it. Hence, be- 
fore applying the formalism of I Lee et ah I (|2006l ). we have 
estimated the strength of Bra emission and the nebular con- 
tinuum in the 4.5 vim IRAC band, and subtracted it from the 
observed flux. 

To estimate the nebular component at 4.5 vim, we first 
calculated the SF R, either from Ha l uminosities whe n avail- 
able, or from L m (|Draine fc Lill2007l ; lKennicuttll 19981 ) if not. 
Th en, assuming the Case B reco mbination coefficients given 
bv lOsterbrock fe Ferlandl |2006l ), we used the observed Ha 
flux (or that inferred from the conversion of SFR to Ha: 
lKennicutt|[l998l ; iLee et al.ll2009l ) to calculate Bra flux. The 
nebular continuum in the 4.5 um IRAC band was inferred 
by interpolating to IRAC wavelengths the volum e emission 
coefficients given bv lOsterbrock fc Ferlandl (|2006l ). 

At low metallicities, T c tends to be rather high, 
£ 10000 K and can even exceed 20000 K (e.g jGuseva et al.1 
l2006l : llzotov et al.ll2007l ). Hence, all calculations were per- 
formed assuming two different electron temperatures, T c , 
10000 K and 20000 K. Because we were unable to uniformly 
derive T c for all the galaxies in our dwarf samples, we sim- 
ply took the average of these two values; this results in an 
intrinsic uncertainty of ~10% for the gas contribution. 

Once the average predicted nebular emission, both Bra 
and continuum, was subtracted from the observ ed 4.5 urn 
photo metry, we applied the method developed by I Lee et al.l 
(2006). This involves a M/L ratio based on the B — K color, 
and a correction to the 4.5 vim band to approximate K, as 
described in Sect. 12. ll 

Because hot-dust emission in low-metallicit y starbursts 
can also be significant at or longward of 2 vim (|Hunt et al.l 
|2001| J 2002), we also attempted to infer the contamination at 
IRAC wavelengths from hot dust. This was possible when 
3.6 vim fluxes were availablfl After subtracting the nebu- 
lar emission from the 3.6 urn fluxes, we then estimated the 
stellar emis si on at 4.5 vim using the formulation given by 
iHelou et al. I (|2004h . i.e., 0.596 times the 3.6 vim emission. 
The excess 4.5 urn emission, after correcting for nebular con- 
tamination and subtracting the stellar component, should be 
largely due to hot dust. 

The amount of contamination in the nearby dwarf 
galaxies, both from nebular emission and from hot dust, 
can be significant. In some sources (notably HS 0822+3542, 
IZwl8, IIZw40, Mrk209, Mrk 1315, SHOC391, Tol65), the 
contamination at 4.5 u.m from Bra alone can exceed 10%. 
Nebular contamination at 4.5 u.m can be as high 25-30%, al- 
though the median value among the dwarf samples is ~10%. 
The hot dust emission at 4.5 vim can also be a relatively large 
fraction of the total: in NGC5253, SBS 0335-052, Haroll, 
Pox 4, the hot dust contamination of the IRAC 4.5 vim band 
is >40%. 

There are 7 BCDs studied by iFumagalli et al.1 (|2010l ) 
that have stellar masses derived from combinations of op- 
tical and near-infrared colors that are in common with our 
sample. In general, the masses agree to with ~30% (0.1 dex). 



This is true only for the LVL and Engelbracht samples. 
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2.7 Caveats 



In the end, our sam ple consists of 16 dlrrs (|Lee et al.| [2006). 
113 LVL galaxies dkennicutt et al.l | 200S ) 89 local star- 

20081 ; [Fumagalli et al] 
Green Peas a nd LCGs 
Jzotov et all 1201 ill, a nd 43 LBGs a t z ~ 1 - 3 dErb et all 
20061; IShaplev et al l I200I , l2005al ; iMaiolino et all l2008i 



burst s and BCDs (lEngelbracht et al.l 
l2Q10D, 



2010; Hunt et al. 



803 



Mannucci et al. 20091 ) . This combined sample of ~1070 



galaxies has stellar masses and metallicities derived with 
heterogeneous methods; these could affect our analysis of 
scaling relations, but as we argue below, do not. 

The z = samples have masses based on the 4.5 pm 
luminosity which minimizes variations in the M/L ratio 
Oun fc Tmll2008h . We have also removed nebular and re- 
combination line emission and applied the B — K co l or cor - 
rection to the M/L ratio as advocated by iLee et all i|2006h . 
Hence, these masses should be sufficiently reliable to com- 
pare with those derived from the more sophisticated SED 
fitting methods applied to the higher-redshift samples. It is 
also true that our sample spans 10 5 in stellar masses, so any 
systematics should be overcome by sheer dynamic range. 

As described in Sect.rj^] for the LVL/11HUGS galaxies 
we used the largest of the UV or Ha-derived SFRs. Relative 
to SFRs inferred from the FUV continuum, SFRs calculated 
with Ha can be underestimated at low SFRs. Although the 
differe nces could be due t o the uncertain extinction correc- 
tions (|Boselli et al.l l2009h , they could also result from ef- 
fects in stochasticity in the number of massive stars in small 
clusters. In any case, they should be directly comparable to 
Hq-inferred va lues for higher SFRs in other samples (e.g., 
ILee et alllioogh . and should not affect our main conclusions 
about metal-poor starbursts. 

The metallicities are potentially more problematic since 
most of the local galaxies and all of the LCGs/Green Peas 
have O/H derived from the direct electronic temperature 
method. On the other hand, the high-z samples rely on 
various strong-line calibrations, so there may be some dis- 
crepancy between the two sets of abund ances. Neverthe- 
less, according to iKewlev fc Ellison! (|2008l ). the largest de- 
viations (up to ~0.7dex) occur at the highest metallicities 
and for the photoionization technique. Typical discrepan- 
cies between the T e and the strong-line methods should be 
<0.15dex, at least for the relative comparison. Since our 
sample is dominated by sub-solar metallicities, and because 
the range in metallicities is so large (~ x 100) , we do not 
expect any problems with the metallicity determinations to 
exert a strong influence on the MZR for our sample. 



3 LOW-METALLICITY STARBURSTS 

Figures [T] [2] and [3] show the trends of SFR, oxygen abun- 
dance, and stellar mass for all samples. Here we focus on 
the outliers from the general relations and identify a "low- 
metallicity starburst" . 

Generally SFR increases with metallicity, although the 
correlation is far from perfect. Inspection of Fig. [T] which 
plots SFR against 12+ log(0/H), shows that in two of the lo- 
cal dwarf samples (LVL, dlrr) there is a clear increase of SFR 
with metallicity. The trend gets rather noisy for the BCD 
samples, and disappears completely for the LCG+Green 
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Figure 1. SFR plotted against nebular oxygen abundance. The 
LVL (11HUGS) galaxies are shown as small open or filled cir- 
cles (filled when SFR^0.6 Mq, yr —1 ), with differen t colors corre- 
sponding to Hubble type T as in ILee et al. I Hool): T ^ 8 cyan, 
5 sg T < 8 blue, 3 sC T < 5 magenta, and T < 3 red. Red X cor- 
responds to the dlrr sample; BCDs are shown as (blue) squares 
(Fumagalli+Hunt samples) and (purple) triangles (Engelbracht); 
Green Peas are given by small (green) filled circles; LCGs by 
+; LBGs at z ~ 1 (Shapley), 2 ~ 2 (Erb, Shapley) are shown 
as 6-pronged asterisks, and at z ~ 3 (Maiolino, Mannucci) as 
8-pronged asterisks. 6-sided open stars show the two dwarf "pro- 
totypes", SBS 0335— 052 and IZwl8. Solid symbols show those 
BCDs with SFR^0.6 Mq yr _1 . The region in the upper left quad- 
rant delineated by dotted lines corresponds to what we loosely de- 
fine as a "low-metallicity starburst" (see text): SFR^0.6 Mq yr _1 
and 12+log(O/H)^8.0. 



Peas. Nevertheless, the statistical significance of the corre- 
lation of SFR with O/H over all the samples (1070 data 
points) is quite high, > 10 a. In fact, two regions tend to be 
underpopulated: the lower right quadrant with high abun- 
dances and low SFRs, and the upper left quadrant with low 
O/H and high SFRs. The galaxies that occupy this latter 
region are what we will call "low-metallicity starbursts". 

Figure [2] shows the "Star Formation Main Sequence" or 
SFMS. In addition to the data, three trends are shown as 
black curves extending only as far as the parameter space in 
the orig inal samples they w ere derived from. The SFMS de- 
rived bv lSalim et all (|2007h for the Local Universe at z ~ 0.1 
(both the linear relation and the "Schechter-li ke" curve) ap- 
pears at lower sSFRs; the one proposed by iNoeske et al.l 
l|2007h for galaxies at redshifts between 0.2 and 0.7 is roughly 
parallel to the previ ous one but with a larger offset; and 
the SFMS derived bv lBauer etail (|201ll ) for a mass-selected 
sample between redshifts 1.5 and 3.0 has an even greater off- 
set. The figure clearly shows that most of the AMAZE and 
LSD samples at z ~ 3 follow the Bauer et al. relation fairly 
well. However, the high-z trend is also reflected in the LCGs 
and the Green Peas at 2 ~ 0.3, and several of the BCDs at 
z ~ 0. Most of the BCDs follow the Local Universe SFMS by 
ISalim et all (|2007l ). but some galaxies at low stellar masses 
have excess SFR relative to the SFMS. At stellar masses 
M star ~ 10 8 M , the SFR of these outliers is > 0.6 Mq yr" 1 . 
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Figure 2. sSFR plotted against stellar mass, M a tar (star- 
formation "Main Sequence"). Symbols are as in Fig. [T] The 
black curve and line correspond to the Local-Universe (z = 0.1, 
ISalim et al J2007h ; the curve reproduces the best-fit Schecter func- 
tion, and the line the best fit to "pure" (no AGN) star-forming 
galaxies. The black line with the higher nor malization gives th e 
intermediate-redshift trend (0.2 < z < 0.7. iNoeske et al.|[2007h ; 
and the bordeaux curves with roun d data points ind icate the 
high-redshift relations (1.5 < z < 3.0. lBauer et al.ll20lih . 
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Figure 3. Nebular oxygen abundance, 12+ log(0/H), plotted 
against stellar mass, M s t a r (MZR, FMR). Symbols are as in 
Fig. [T] The blac k dotted line gives the linear relation found by 
iLee et alj (2006), and t he curves correspond to the adjustments 
for SFR as advocated bv lMann ucci et al.| l|201ll). SFRs consis tent 
with the ranges considered by |]vEmnuTciet*lltL | j2010l, 1201 ll) arc 
shown (0.1 Mq yr -1 , dashed line to the left; 1 Mq yr — 1 , solid line; 
10 Mq yr -1 , dotted line to the right). 



This motivates our choice of SFR = O.6M0yr as a SFR 
threshold for starbursts at low metallicity. Such a low-mass 
galaxy would have a n sSFR of ~ 6 x 10 ~ 9 yr" 1 , well outside 
the range probed bv lSalim et alj (2007), but roughly consis- 
tent with an extrapolation to lower masses of the "starburst 
sequence" they define. 

The MZR for our data is shown in Figure [3] The 
FMR proposed by M annucci et alj l|201Cf ) is also shown 
as a set of black curves; these curves are the predictions, 
based on the SDSS, of how SFR alters the MZR. The 
dashed curve corresponds to SFR = 0.1 Mq yr -1 , i.e., where 
a galaxy with this SFR would lie relative to the MZR. The 
solid curve gives a correction for SFR= 1 Mq yr -1 , and the 
dotted one SFR = 10 Mq yr " 1 ; the curves span rou ghly the 
ranges in parameters of the lMannucci et alj l|2010t) sample. 
The FMR was s ubseq uently extended to lower masses by 
iMannucci et al.1 (|201ll ) (shown here as the linear extensions 
toward lower masses slightly detached from the m ain curved 
FMR) . The linear MZR for low masses found by ILee et al.l 
(2006) is also shown, but th e slope with M B ta r is flatter 
(~0.3) than the slope found bv lMannucci et al.l (|201ll . ~0.5). 

Most of the AMAZE and LSD galaxies at z ~ 3 follow 
neither the MZR nor the FMR , and this was already no- 
ticed bv lMannucci et all l|2010l ). However, some of the low- 
metallicity dwarf galaxies with SFR>0.6 Mq yr -1 are also 
outliers of both the MZR and the FMR, although galax- 
ies with lower SFR follow the MZR and/or FMR quite 
well. The degree of deviation from the main trends increases 
for metal abundances at or below 12+ log(O/H)~8.0. Some 
metal-poor BCDs and high-z LBGs in the AMAZE+LSD 
samples deviate from the MZR and FMR by factors of 100 
or more in stellar mass. Hence, we adopt this abundance 
12+ log(0/H) = 8.0 as a "fiducial location" in the MZR be- 
low which we can distinguish potential outliers; here we are 
interested in low-metallicity starbursts, and consequently 
concentrate on this metal-poor abundance regime and on 
high SFRs. The region delineated in Fig. [1] is occupied by 
these metal-poor starbursts, independently of their redshift. 
If a galaxy has 12+ log(O/H)^8.0 and SFR>0.6 Mq yr" 1 , 
we will call it a "low-metallicity starburst" (LMS). 

3.1 Crude demographics 

The only galaxy sample in our database which can be con- 
sidered complete in any sense is the 11HUGS/LVL sample; 
it is approximately distance (thus volume) limited, and a 
fair representation of the nearby galaxy population. Of 113 
galaxies in the LVL sample, only 2 (NGC 784 and NGC 4656) 
would be considered LMSs, a fraction of ~ 2%. The fraction 
of low-metallicity starbursts in the BCD samples is 5 times 
higher, ~11% (10 LMSs), but these galaxies have already 
been selected from the general galaxy population through 
the strength of their emission lines. 

Among the LCGs and Green Peas, the fraction is higher 
still, but, again, these galaxies have been selected for strong 
emission lines; 176 galaxies out of 803 LCGs would be clas- 
sified as a LMS, roughly 22%. This is partly a result of a 
selection effect because of the limit on H/3 luminosity, and 
thus SFR, imposed by llzotov et al] l|201lf ); all LCGs have 
SFRs greater than our fiducial cutoff by definition. At z ~ 3, 
there are 3 LMSs in the LSD sample of 9 and 2 LMSs of 9 
AMAZE galaxies, or ~27%; these are small-number statis- 
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tics, but roughly consistent with the LCGs at a redshift 10 
times smaller. 

Outside the LVL and dlrr samples, virtually all the 
galaxies are selected on the basis of high equivalent widths 
in various emission lines. These samples have high LMS frac- 
tions, perhaps more similar to the LBGs at z ~ 3 than either 
the LVL or dlrr samples. This means that the differences 
in properties found for high-z populations may result from 
differences in selection criteria, relative to the local galaxy 
populations. We may be observing locally a somewhat rare 
class of low-metallicity starbursts that are however similar 
to some high-z galaxy populations such as LBGs. In the 
next section, we put this proposition on a firmer footing by 
showing the existence of a Fundamental Plane of O/H, SFR, 
and M s t ar , which apparently does not vary with redshift. 



4 THE FUNDAMENTAL PLANE OF 

METALLICITY, SFR, AND STELLAR MASS 

Unlike the SE > SS s a mples stud i ed b y other groups 



llTremonti et al.1 |2004 ISalim et"al] 120071 ; IMannucci et all 
|2010| ; lYates et alj|2012h , our sample has been designed to 
probe low stellar masses, low metallicities, and (relatively) 
high SFRs. This makes it uniquely suited to exploring a 
much expanded parameter space, and the resulting scaling 
relations of metallicity, M a t ar , and SFR. At high masses and 
metallicities, both the MZR and SFMS inflect and flatten. 
Howev er, for M Btar ^ 3 x 10 10 Mp), roughly the "turn-over 
mass" |Tremonti et aij|2004 IWvder et alj|2007l 1. the trends 
are very close to linear, independently of redshift. Hence, 
the observationally-defined variables (O/H, SFR, and M star ) 
should define a plane which, given the relatively large scat- 
ters in the SFMS, the MZR, and the SFR-O/H relati on, is 
not optimally viewed (see also lLara- Lopez et al.|[20ich . 

We can implicitly test the assumption that the relations 
among the observationally-defined variables are invariant, 
regardless of redshift. Hereafter, for brevity, we will call the 
"observationally-defined variables" observables, even though 
they are derive d from observatio n s rath er than being directly 
observed. Like IMannucci et al.l I2010T ). we propose that if 
the MZR and SFMS are observed to "evolve" or change 
with redshift, then this evolution or alteration must be due 
to different selection criteria for diffe rent galaxy popu l ation s 
at different lookback times. However. IMannucci et all (|2010T ) 
defined a "surface" by associating SFR with stellar mass, 
and plotting this combination against metallicity; they did 
not examine the mutual correlations among the three pa- 
rameters. Metallicity, mass, and SFR are all correlated (see 
Sect. but it is not clear which parameters are mainly 
responsible for the dispersion. 

We have therefore performed a Principle Component 
Analysis (PCA) on all galaxies with stellar masses ^ 3 x 
10 10 M Q (1019 data points). A PCA diagonalizes the co- 
variance matrix, thus giving the linear combinations of ob- 
servables which define the orientations that minimize the 
covariance. These orientations are the eigenvectors and are, 
by definition, mutually orthogonal. In the 3-space formed by 
the three observables, to form a plane we would expect most 
of the variance to be contained in the first two eigenvectors; 
for the third, perpendicular, eigenvector the variance should 
be very small. 



Our sample is particularly well suited for such an anal- 
ysis; it spans almost two orders of magnitude in metallic- 
ity (12+log(0/H) = 7.1 to ~ 9), a factor of ~ 10 6 in SFR 
(~ 10~ 4 < SFR ^~ lO 2 M yr _1 ), and a factor of ~ 10 5 
in stellar mass (~ 10 6 < Af star <~ 10 n M Q ). Other sam- 
ples previously analyzed to find scaling relations cover much 
smaller parameter ranges: typically less than a decade in 
metallicity (12+ log(0/H)^8.4), a factor of ~200 in SFR 
(~0.04< SFR< 6M Q yr" 1 ), and roughly 2 orders of magni 



tude i n stellar mass (M a t ar > 10 Mq) iTremonti et al 

| 2004l ; ISalim et al] 120071 ; IMannucci et al.l I2010L I Yates et al 
|2012j). This implies that our results should represent the 
extremes of SF processes, possibly closer to what is ob- 
served in the early universe than nearby. To avoid exces- 
sively weighting the LCGs, which effectively dominate our 
sample (~800/1000), we used a bootstrapping technique to 
randomly select a fraction of LCGs and repeat the PCA 
many times on different samples of ~350-1000 objects. 

The PCA shows that the set of three observables truly 
defines a plane; ~97-98% of the total variance is contained in 
the first two eigenvectors. 80-85% (depending on the boot- 
strapping) of the variance is contained in the first eigen- 
vector or Principle Component (PC), PCI; 12-17% of the 
variance is found in the second PC, PC2; and <2-3% in 
PC3. Hence, the 3-space defined by O/H, SFR, and M star is 
degenerate; only two independent parameters are necessary 
to describe the galaxies in these three dimensions. 

The existence of a plane of O/H, SFR, and M star 
for star- forming galaxi es is not, however, a new result. 
lLara- Lopez et ahl ( 2010T I found that a plane could describe 
a sample of ~ 33 000 SDSS galaxies but, rather than per- 
forming a PCA, they fitted linear regressions to pairs of 
observables. They also included a second-order polynomial 
for low-metallicity mass bins, and concluded that a "lin- 
ear combination of the metallicity and SFR as a function 
of stellar mass" could describe their sample. Thus, unlike 
IMannucci et al.1 d201Ch who combined SFR a nd M s t ar to pre- 
dict 12+ log(Q/H), lLara-L6pez etafl (|2010h combined SFR 
with 12+ log(0/H) to predict M star , and find a dispersion a 
of 0.16 for the 33 000 SDSS galaxies in their sample. 

For our sample, Fig. [4] shows three different projections 
of the plane defined by these three "basis" (eigen-) vectors 
that are different linear combinations of the observed vari- 
ables: 

PCI = 0.12a;i +0.75 z 2 + 0.65x3 (1) 
PC2 = -0.3121-0.65x2-0.69 23 (2) 
PC3 = -0.9421 -0.11 x 2 + 0.31 23 (3) 

The PCA first normalizes the variables to zero mean, so that 

21 = 12 + \og{0/H) - (12 + \og{0/H)) 
= 12 + \og(0/H) - 8.063 

22 = \og(SFR) - (\og{SFR)) 

= log(SFR) + 0.594 M yr" 1 

23 = log(Mstar) - (log(Mstar)) 
= log Mstar -8.476 M Q 

The top panels of Fig. [4] show the best edge-on views of the 
plane, and the bottom panel its face-on view. 

The first eigenvector PCI is weighted roughly equally 
between SFR and M star ; its dependence on O/H is relatively 
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Figure 4. Different projections using the 3 PCs found by the PCA. Galaxies are coded as shown in the lower right (empty) panel. 
The top left and right panels show the orthogonal "edge-on" views of the plane; the bottom panel shows the plane face-on. x\ = 
12+log(0/H)-{ 12+log(0/H)), x 2 = log(SFR) — (log(SFR)) (M yr" 1 ), and x 3 = log(Af star ) - <log(A/ star )) (M ). (12+ log(0/H)) = 
8.063; (log(SFR)}= —0.594; (log(M star )} = 8.476. PCI = 0.12 Xl + 0.75 x 2 + 0.65 x 3 ; PC2 = -0.31 Xj - 0.65 x 2 - 0.69 x 3 ; PC3 = 
-0.94a;i - 0.11x2 + 0.31x 3 . 



small. The same is true for PC2 with the relative weights be- 
tween SFR and M sta r reversed, but with a slightly increased 
dependence on O/H. This means that most of the dispersion 
in the 3-space is governed by the spreads in stellar mass and 
SFR, rather than by metallicity. 

The eigenvector with the least variance, PC3, is domi- 
nated by metallicity with a secondary dependence on M s t ar . 
Of the three parameters, metallicity is the least influential on 
the sample variance. Apparently, the plane is defined mainly 
by SFR and stellar mass A/ S t ar and its vertical extent or 
thickness by metallicity. However, this result may stem from 
the breadth and coverage of parameter space by our sam- 
ple. In fact, our bootstrapping experiments suggest that the 
PCA coefficients depend slightly on the composition of the 
sample. The sample spread in SFR is particularly important 
as the SFR coefficients can change by 20-30%. Repeating a 
PCA on different samples would be worthwhile in order to 



gauge the generality of our results and the reliability of the 
coefficients (see below). 

Since PC3 contains only ~2% of the variation among 
the galaxies, it is a potentially powerful tool for establishing 
an optimized view of the parameter space defined by metal- 
licity, stellar mass, and SFR. The tightness of PC3 and its 
dominant dependence on metallicity makes it possible to 
formulate a correlation between O/H, and some linear com- 
bination of the other two variables. The PCA has essentially 
distilled the information given by the MZR and the SFMS 
(and the trend of SFR with O/H), thus providing the com- 
bination of the relations with the lowest dispersion. If we 
take PC3 as shown in Fig. [4] and set it equal to zero, we 
can derive a best-fit relation among O/H, SFR, and M star , 
which is characterized by small variation over a wide range 
of parameter space. This is shown in the left panel in Fig- 
ure[5]where we have plotted 12+log(0/H) vs. the equation 
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-0.12 Log(SFR) + 0.33Log(M,) + 5.20 12 + Log(0/H) - FP 

Figure 5. Left panel: the PC3 component shown as 12+log(0/H) vs. the resulting linear combination of (log) SFR and Mstar obtained 
from the PCA on our sample. Galaxies are coded as in Fig. [4] The FP applied to the FMR-SDSS sample is shown as red contours, 
with the innermost contour corresponding to 60% of the sample of 55 082 galaxies, and increasing to 70%, 80%, 90%, 95%, 99%, with 
99.5% as the outermost contour. The dotted line in the left panel indicates equality. Right panel: the histogram of the residuals from the 
identity (dotted line) relation shown in the left panel; the blue Gaussian and broader distribution corresponds to our sample, and the 
red Gaussian and narrower distribution to the FMR-SDSS galaxies. The FMR-SDSS distribution is well fit by a Gaussian but offset by 
+0.4 dex. In both samples (and both panels), only galaxies with M s tar < 10 10 ' 5 Mq are shown. 



that results from equating PC3 to zero (only galaxies with 
Mstar < 3 x 10 10 M Q are shown). 

12+log(0/tf) = -0.121og(SFR)+0.331og(M star ) + 5.20 (4) 

The equality line in the left panel and the residuals in the 
right panel show that the observed abundance is well esti- 
mated by the FP. The reliability of this estimate is roughly 
constant over the entire parameter space covered by our 
sample; the errors shown in the right panel of Fig. [5] are well 
approximated by a Gaussian with a a = 0.17, or <48-50% 
uncertainty (0.17dex). 

Also shown in Fig. [S] is t he FP applied to the SDSS 
sample used to define the FMR (|Mannucci et al.ll201fj ); only 
galaxies with Mstar ^ 3 x 10 10 Mq are plotted. The disper- 
sion of the SDSS sample around the FP defined h ere is small, 
a = 0.06 dex, the same as lMannucci et"aL I (|2010l ) found rela- 
tive to the FMR. However, the SDSS sample appears to be 
offset from the FP by ~0.4dex, as shown in the right panel 
of Fig. The mean metallicity of the FMR-SDSS sample 
(with 55 082 galaxies having masses below 10 10 - 5 M Q ) is 8.94 
(±0.14), roughly the limit of t he maximum oxygen abun- 
dance in spiral galaxies (8.95: iPilvugin et alj [20071 ). Such 
metallicities are also well into the metallicity range known 
to have large d iscrepancies ( 2i 0.4 dex ) relative to direct ca l- 
ibrations (e.g-. lKewlev fc Ellisonll2008l ; |Pilvugin et al.ll2010l ). 
This is not a problem peculiar to the SDSS sample, since the 
metal-rich galaxies we examine here fall in a similar metal- 
licity region (see F ig. [5]). 

As stated by fMoustaka s et all |2010l ). "the factor of 
5 absolute uncertainty in the nebular abundance scale 



poses one of the most important outstanding problems in 
observational astrophysics". There is some evidence that 
calibr ations based on theo re tical strong-line calibrations 
(e.g., ITremonti et all 12004 iKobulnickv fc Kewlevl |2004 
iNagao et al.ll2006l ) overestimate the n ebular metallicity by 
0.3-0.5 dex (e.g.. IBresolin et al.l l2009), but direct methods 
may underestimate it by 0.2— 0.3 dex (e.g.. [Przvbilla et al.l 
2008). Such discrepancies may be due to temperature mix- 
ing i n the Hi I regio ns in metal-rich giant spiral galaxies 
(e.g., |Peimbert|[l967l), but are much less common in metal- 
poor dwarfs (e.g.. |Pilvugin et al.ll201ol ). In any case, different 
metallicity calibrations make comparing samples over wide 
ranges of abundances quite problematic, and are ultimately 
a major obstacle to probing scaling relations over a vast pa- 
rameter space. 

The FP fits the FMR-SDSS sample as well as the FMR 
itself (<r =0.06 dex), although with a +0.4 dex offset which 
could be due to proble ms with the metallici ty calibration 
at high abundances (e.g.. IBresolin et al.ll2009l ). The FP dis- 
persion o f 0.17 dex for o u r sam ple(s) is higher than that 
found by ITremonti et~al] (|2004l ) for the MZR defined by 
53 000 galaxies from the SDSS (0.1 dex), and also higher 
than t he FMR (0.06 d e x) fou nd for the similar SDSS sam- 
ple bv lMannucci et all (|201ll ). However, both SDSS samples 
span limited ranges in O/H, SFR, and stellar mass relative 
to the parameter space covered by our sample. Indeed, the 
mean and standard deviation of (log) stellar mass for the 
FMR-SDSS sample (with Af star < 10 10 5 M ) is 10.07+0.31 
(the comparable mean, standard deviation in our sample is 
8.48+1.10). Thus, a lower dispersion in the fit is not surpris- 
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ing despite the many more galaxies in those samples. The 
value of the FP dispersion for the ~1000 galaxies studied 
here is only slightly higher than t hat found for th e MZR of 
25 nearby dwarf galaxies (0.12 dex. lLee et al. | 2~006 | ). another 
sampl e do minated by low - mass galaxies. iTremonti et al] 
(2004) and lMannucci et al.l <|2010h found that dispersion in 
the MZR increases with decreasing mass; hence given the 
parameter space covered by our sample, the a of 0.17 dex is 
reasonable. 

It is also true that our sample is rather heterogeneous, in 
terms of the methods with which metallicity and stellar mass 
are calculated. Unfortunately this is of necessity since we 
are comparing properties of local and high-redshift galaxy 
populations, over a wide range of metallicity. Thus, it is 
possible that intrinsic systematics could be responsible for 
a fraction of the dispersion observed for our sample. 

4.1 PCA of the SDSS sample 

Because our original aim was to preferentially select metal- 
poor starbursts, our combined sample could be biased by 
metallicity. In such a case, the FP found here would not 
be applicable to general galaxy populations. To test this 
po ssibility, we performed a PCA of the SDSS sample used 
by iMannucci et alj |2010f) , considering only those galaxies 
(55 082) with M st ar< 3 x 1O 1O M , the same limit as im- 
posed for our sample. This is the same FMR-SDSS sample 
discussion in the previous section. If the coefficents and dis- 
persions are similar to the ones obtained from the PCA on 
our sample, we could conclude that our sample is not unduly 
biased, and that the results can be generalized to broader 
galaxy populations. 

The eigenvectors resulting from the PCA on the FMR- 
SDSS are found to be very similar to those we obtained 
for our sample. Specifically, 98% of the variation among the 
parameters is contained in the first two eigenvectors, PCI 
and PC2, and they are similarly dominated by M sta r and 
SFR. Moreover, PC3, dominated by O/H, comprises only 
2.3% of the variance. Thus, following the same reasoning 
as before, we have set PC3 equal to zero and solved for 
12+log(0/H). From the PCA on the 55 082 galaxies in the 
mass- limited FMR-SDSS, we find: 

12+log(0/tf) = -0.121og(SFR)+0.381og(M Bta r)+5.06 (5) 

Comparison with Eq.|3] shows that the coefficient multiply- 
ing (log) SFR is the same as for our sample, but there is 
an offset of —0.14 in the constant. The M sta r coefficient is 
also slightly steeper than for our sample (0.38 vs. 0.33), al- 
though not as steep as the 0.51 coefficient multiplying M ata ,r 
in the extension of the FMR to low mass bv lMannucci et ahl 
(|201lT ) . The difference in the FP constants (-0.14) is not 
as large as the offset between the two samples (—0.4, see 
Fig. [5]), but the larger coefficient multiplying (log) M sta r 
compensates the difference. Fig. [S]shows both samples plot- 
ted with the FP defined by the FMR-SDSS galaxies, and 
the relative dispersions. 

If we apply this new (SDSS) FP to our sample, we find 
the same spreads for both samples as with the original FP 
(0.17dex, 0.06 dex, see Fig.[6|. However, our sample requires 
an offset of —0.3 dex, slightly smaller (and of opposite sign 
as expected) than the SDSS offset relative to the original 
FP (+0.4 dex). Hence, the two formulations are probably 



equally good representations of galaxies up to z ^ 3, across 
a wide range of stellar mass (10 6 <M star < 3 x 1O 1O M ), 
metallicity, and SFR (10 -4 <SFR< ^Mgyr" 1 ). Never- 
theless, as discussed above, differences in metallicity calibra- 
tions make it difficult to accurately assess scaling relations 
across wide ranges of abundances. 



4.2 Comparison with the FMR 

As mentioned in the Introduction, M annucci et all (|2010h 
devised the Fundamental Metallicity Relation which sig- 
nificantly reduces the scatter in the MZR of SDSS- 
selected galaxies by in troducing a dependence on the SFR. 
IMannucci et al.l ()201lh extended the FMR to lower masses 
(^0.32 < 9-5) such that: 

12 + log(0/H) = 0.51(^0.32 - 10) + 8.93 

= -0.16, log(SFR) + 0.51 log(M star ) + 3.83 (6) 



where ^0.32 = log(M star )-0.32 log(SFR). 

We have applied the FMR to our sample (considering 
also the non-linear extension to high masse^ 4 ]) and find that 
the residuals have mean values of >0.32 dex, almost twice 
the value we find here for either FP. The FMR applied to 
our sample and to the FMR-SDSS sample is shown in Fig. (7) 
The residuals for our sample are asymmetrically skewed to 
negative values, and fitting them to a Gaussian gives a = 
0.27 dex. Hence, the FMR is a poorer fit to our sample than 
the FP. 

The reason is that th e exten sion to lower masses pro- 
posed by IMannucci et "all |201ll 'l. with a slope of 0.51, is 
too steep to fit the galaxies in our sample (see also Fig. [T]). 
From the PCA, we find a meta llicity-mass dependence of 
0.33, similar to iLee et alj (|200fJ ) who found for their low- 
mass dlrr sample a linear dependence (in log-log) of 0.30 
between 12+log(0/H) and M star (a =0.12 dex). This shal- 
lower dependence of metallicity on mass may be peculiar to 
samples with many low-mass galaxies, although we find sim- 
ilar values even when restricting the fits to only the highest 
masses, and when fitting the FMR-SDSS sample. 

Nevertheless, the relative dependence between M atar 
and SFR found by o ur PCA is similar to that given by 
IMannucci et all l^Old ). With Mo. 32, they assumed a unit de- 
pendence of O/H on (log)M s t ar , and found that the value 
of 0.32 multiplying (log)SFR minimizes the residuals of 
the MZR relative to u . The extension to lower masses by 
IMannucci et al.1 |201ll ) has a flatter dependence of O/H on 
Mo. 32, 0.51 (see above), but the same value of 0.32 relat- 
ing (log) SFR to M star . For the FP, we find a shallower de- 
pendence of O/H on M s tar (0.33), and a SFR coefficient of 
—0.12; multip lying 0.33 (our mass d ependence) by the 0.32 
SFR factor of IMannucci et al.l (|2010h would give a SFR de- 
pendence of —0.11, very similar to the value of —0.12 given 
by the PCA. 



4 12+log(O/H) = 8.90+0.37m-0.14s-0.19m 2 +0.12ms-0.54s 2 
for no.32 ^ 9-5, where m = log(M B tar) — 10 an d s = log(SFR) in 



solar units. See Eq. 2 in IMannucci et all i201lf ). 
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Figure 6. Left panel: the PC3 component shown as 12+log(0/H) vs. the resulting linear combination of (log) SFR and M s ta.v obtained 
from the PCA on the FMR-SDSS sample. Galaxies are coded as in Fig. [4] As in Fig. [5] the FMR-SDSS sample is shown as red contours, 
with the innermost contour corresponding to 60% of the sample of 55 082 galaxies, and increasing to 70%, 80%, 90%, 95%, 99%, with 
99.5% as the outermost contour. The dotted line in the left panel indicates equality. Right panel: the histogram of the residuals from the 
identity (dotted line) relation shown in the left panel; the blue Gaussian and broader distribution corresponds to our sample, and the red 
Gaussian and narrower distribution to the FMR-SDSS galaxies. The distribution of residuals for our samples is well fit by a Gaussian 
but offset by — 0.3dex. In both samples (and both panels), only galaxies with A/star < 10 10 ' 5 Mq are shown. 




FMR 12 + Log(0/H) - FMR 

Figure 7. Left panel: nebular oxygen abundance, 12+ log(Q /H), plotted against t he FMR abscissas, as described in the text (there are 
two separate branches, according to the value of /io.32i see iMannucci ct al.|[201lf) . Symbols and contours are as in Fig. [5] the dotted 
line indicates equality. Right panel: the histogram of the residuals from the identity (dotted line) relation shown in the left panel; the 
blue Gaussian and broader distribution corresponds to our sample, and the red Gaussian and narrower distribution to the FMR-SDSS 
galaxies. Our sample is offset by — 0.18dex, with a slightly asymmetric distribution. In both samples (and both panels), only galaxies 
With Mstar < 10 10 - 5 Mq are shown; this mass limit is clearly evident in the left panel. 
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5 DISCUSSION AND CONCLUSIONS 

From a sample of ~ 1100 galaxies spanning a wide range 
of O/H, SFR, and Mstar, we have identified a class of 
low-metallicity starbursts which deviate from the SFMS 
and the MZR in the same way, independently of redshift. 
When these scaling relations are considered together, the 
sample as a whole (although excluding the galaxies with 
M s tar> 3x 10 10 Mq) defines a Fundamental Plane whose ori- 
entation is defined primarily by SFR and Mstar, and whose 
thickness is governed by O/H. 

The two-dimensional nature of the plane in the 3-space 
defined by the observables implies that on ly two parameters 
are n ecessary to describe a galaxy (see also lLara- Lopez et al.l 
2010). The thinness of the plane when viewed edge-on en- 
ables the formulation of a relation which combines the SFMS 
and the MZR in an optimal way. Through application of a 
PC A, we have established a dependence of O/H on SFR 
and M s t a r which has a deviation over the whole sample (for 
M st ar< 3 x 10 10 M©) of 0.17dex, or ~48%. 

The FP from the PCA applied to the S PSS sample 
used to define the FMR (Ma nnucci et alj|20ich is similar to 
that obtained for our sample (see Eqs. 0] [SJ . Applying ei- 
ther FP to our sample and to the FMR-SDSS shows an offset 
in metallicity which distinguishes the two, but the disper- 
sions are similarly low (0.17 dex and 0.06 dex for our and 
the FMR-SDSS samples, respectively). However, the FMR 
is not the best description of our sample, as it gives a higher 
dispersion (0.27 dex) than either FP. 

Perhaps the most intriguing aspect of the FPs given 
by the PCA is that they hold for all the redshifts we have 
examined, z < 3. In our sample, the LBGs at all redshifts 
follow the "Fundamental Plane" relation which combines 
the MZR and the SFMS. There are no high-redshift outliers. 
Moreover, the low-metallicity starbursts at all redshifts also 
follow the relation. On the basis of O/H, SFR, and stellar 
mass alone, apparently there is no distinction between the 
behavior of galaxies in the Local Universe and those at high 
redshift (at least up to < 3). 

The statistical analysis does not clarify which relation 
or which observable is the most important; it merely in- 
dicates that there is a plane that describes the family of 
galaxies in the 3-space of O/H, SFR, and M sta r- However, 
the small amount of variance contained in the eigenvector 
dominated by metallicity (PC3) suggests that metallicity is 
the least important of the three parameters in driving the 
variance of our sample. The physics of why or how these 
parameters conspire to form a plane in dependent of redshif t 
will be explored in a companion paper (IMagrini et al.ll2012l ). 
There we show that the deviations from the main scaling re- 
lations are due to a different mode of star formation, a star- 
burst or "active" mode which evolves over relatively short 
timescales. The main scaling relations are instead defined by 
a more quiescent, "passive" mode of star formation which 
evolves more slowly. 

The existence of a Fundamental Plane that spans red- 
shifts up to z > 3 implies that the SFMS and MZR do 
not really change with lookback time, but rather that the 
galaxy populations that define them change with redshift. 
Starbursts become increasingly common at high redshift, 
but they also exist locally and at low metallicity. Thus, the 
"evolution" in the scaling relations is a result of selecting 



those galaxies which are most common at a particular red- 
shift, and may not reflect a true change in the physics of 
how galaxies evolve. 
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